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In the present contribution we review basic mathematical results for three physical sys- 
tems involving self-organising solid or liquid films at solid surfaces. The films may undergo 
a structuring process by dewetting, evaporation/condensation or epitaxial growth, respec- 
tively. We highlight similarities and differences of the three systems based on the observa- 
tion that in certain limits all of them may be described using models of similar form, i.e., 
time evolution equations for the film thickness profile. Those equations represent gradient 
dynamics characterized by mobility functions and an underlying energy functional. 

Two basic steps of mathematical analysis are used to compare the different system. First, 
we discuss the linear stability of homogeneous steady states, i.e., flat films; and second the 
systematics of non-trivial steady states, i.e., drop/hole states for dewetting films and quan- 
tum dot states in epitaxial growth, respectively. Our aim is to illustrate that the underlying 
solution structure might be very complex as in the case of epitaxial growth but can be better 
understood when comparing to the much simpler results for the dewetting liquid film. We 
furthermore show that the numerical continuation techniques employed can shed some light 
on this structure in a more convenient way than time-stepping methods. 

Finally we discuss that the usage of the employed general formulation does not only relate 
seemingly not related physical systems mathematically, but does as well allow to discuss 
model extensions in a more unified way. 

The paper is published in: /. Phys.-Cond. Mat. 22, 084019 (2010) and 



can be obtained at http://dx.doi.Org/10.1088/0953-8984/22/8/084019 



* homepage: http://www.uwethiele.de 
^Electronic address: u.thiele@lboro.ac.uk 



2 



I. INTRODUCTION 

Structure formation at interfaces and surfaces occurs widely in our natural and technological 
environment. The spectrum of related phenomena ranges from growing dendrites in solidification 
or crystallization to budding membranes in the biological cell. Technological processes based 
on or affected by interfacial structuring processes involve, for instance, the sputtering of solid 
surfaces (that may roughen), the usage of instabilities in the epitaxial growth of nano- or quantum- 
dots, the structuring of homogeneous liquid or elastic coating layers, the deposition of structured 
nano-particle assemblies employing instabilities, and heat-exchanger technology based on transfer 
enhancement by surface waves on falling liquid films. A selection is discussed in [TJ. 

Part of the mentioned structuring processes can be modeled as an evolution in time of a surface 
profile. Models are normally based either on a stochastic microscopic discrete or a deterministic 
mesoscopic or macroscopic continuum approach. For an overview of methods for the modelling 
of (solid) nanostructures see fl2j. Here, we focus on the continuum approach. There exists an 
important subset of systems that evolve towards an equilibrium state corresponding to an energetic 
minimum, i.e., normally these are relaxational systems without any external forcing. 

A continuum description of relaxational systems can often be brought into the form of time 
evolution equation(s) for one or several conserved or non-conserved order parameter fields 0(x, t) 
(cf. [3]). A non-conserved field might still follow in part a conserved dynamics. The dynamics is 
governed by the underlying energy functional F[<f)\. The simplest form for a time evolution of a 
purely dissipative system without any inertial hamiltonian dynamics corresponds to the gradient 
dynamics 

d t 6 = V- 
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with the M c (<p) > and M nc (0) > being the mobility functions for the conserved and non- 
conserved part of the dynamics, respectively. Here and in the following d t and d x denote partial 
derivatives w.r.t. time and space, respectively. 

A typical example is the Cahn-Hilliard equation describing the demixing of a binary mixture, 
i.e., a purely 'conserved dynamics' (M nc = 0) ||3]-0. Another example is the Allen-Cahn equation 
describing, for instance, the dynamics of the Ising model in the continuum limit [3]. Multiplying 
Eq. ([I]) by 5F/5(f> and integrating one obtains after a partial integration 

<*W- UJv'lXiV- lMj S lXtV<V (2) 
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confirming F[<f>] to be a Lyapunov functional. 

In the context of evolving surfaces or interfaces, such an equation appears in various contexts. 
We will discus here (i) film thickness equations for films of non-volatile and volatile liquids on 
solid substrates and (ii) surface profile equations for epitaxial growth. 



liquid film on a solid substrate under the influence of capillarity and wettability. In this case, the 
function 0(x, t) represents the film thickness profile and the functional F[<j)] is given by 



where 7 is the liquid-gas surface tension and f(<p) is a local free energy, related to the disjoining 
pressure 11(0) by II = —d^f((f>). The term //</> represents an overall energy bias towards the liquid 
or the gaseous state. It is the sole responsible for evaporation/condensation of flat 'bulk' films, 
i.e., films that are thick as compared to the range of the disjoining pressure. \i corresponds to a 
chemical potential. For details and specific choices for the disjoining pressure see, e.g., Refs. [|6]- 



For Poiseuille flow in the film without slip at the substrate the mobility for the conserved part 
is M c = (j) 3 /3r] where i] is the dynamic viscosity. Several slip regimes might be accounted for by 
different choices for M c ((/>) lfT2ll . The mobility function for the non-conserved part is normally 
assumed to be a constant (see |[T3~l0 . Note that there exists an ongoing discussion regarding the 
form of the non-conserved part of the dynamics (cf. e.g. 1173141810 . 

Eq. ([[]) with ([3]) is extensively studied in the conserved case, i.e., for non-volatile films (M nc = 
0). It can easily be derived from the Navier-Stokes equations and appropriate boundary conditions 
at the substrate and the free surface employing a long-wave or lubrication approximation [|9l ITTTl . 
Depending on the particular physical situation studied many different forms for the local energy 
function /(</>) are encountered. Beside 'proper' disjoining pressures that model effective molecular 
interactions between film and substrate (wettability) j6l Ul OH [I9l - l24ll the equations may as well 
incorporate other pressures modelling, e.g., the influence of an electric field on a film of dielectric 
liquid in a capacitor ll25l - |29l or films on homogeneously heated substrates that form structures 
due to a long-wave Marangoni instability lTT4ll30] - [34ll . The latter is especially interesting because 
it represents a system that is kept permanently out of equilibrium but is nevertheless described 
by a gradient dynamics. Note that the situation is slightly different in a closed two-layer system 
ll3~5ll that can in two dimensions be described by Eq. ([[]) with an appropriate F[(p) but not in three 



Eq. ([[]) might describe the evolution of the surface profile of an evaporating or condensing thin 




(3) 
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dimensions. 

The above mentioned Cahn-Hilliard equation describing the conserved dynamics of demixing 
of a binary mixture corresponds to Eq. ([T]) and Eq. ^ with a constant M c , M DC = and /(</>) being 
a symmetric double well potential. In consequence, many results obtained for the decomposition 
of a binary mixture have a counterpart in the dewetting of thin films and vice versa. The analogy 
was first noted by Mitlin [8J resulting in the notion of 'spinodal dewetting'. Note that there exist 
other choices for /(</>) and M c in the Cahn-Hilliard equation Il36ll37ll . 

Equations of similar form may as well model the epitaxial evolution of surfaces of crystalline 
solids [|2l 13814401 . We illustrate this employing one of the local models for Stransk-Krastanov 
growth found in the literature - namely a simplified 'glued wetting-layer model' (for details and 
derivation see refs. BUI eqs. 18-20] and ETJ). The model assumes an isotropic wetting energy 
of the epitaxial film on the solid substrate, that is added to the (anisotropic) surface energy. Elas- 
tic stresses act through a destabilizing surface stiffness term. It is furthermore assumed that a 
given amount of material is deposited on the surface that then re-arranges in a process of self- 
organisation that might lead to the creation of nano- or quantum-dots, i.e., localized surface struc- 
tures on the nanometre lengthscale. In the case of high-symmetry orientations of a crystal with 
cubic symmetry the evolution of the surface profile 0(x, t) is described by Eq. ([I]) when using 
small-slope approximation. As the amount of material is fixed only the conserved part contributes, 
i.e., M c > and M nc = 0. The model in [|40l employs a constant mobility M c (corresponding to 
a constant surface diffusion coefficient). Non-constant mobilities might as well be used. Note that 
a fourth-order kinematic term is omitted in the evolution equation as it can lead to artifacts if the 
slope of the interface is large (inside the small slope approximation, for details see pT|). The free 
energy functional is 



|(V0) 2 + ^(A0) 2 + ^(V0) 4 + /(0)- 



dV (4) 



where a > is the destabilizing surface stiffness resulting from elastic stresses, v > represents 
the energetic cost of corners and edges, and a quantifies the slope-dependent anisotropic surface 
energy (note that we fixed the b of [|40l as b = a/3 to simplify the equation for the present purpose 
of comparison). The local free energy /(</>) = J W Q ((f))d(f) results from the wetting interaction. 
Thereby W = —w((j)/5Y aw exp(— <p/5) where 8 is a characteristic wetting length, the positive 
w and a w characterize the strength and singularity of the underlying interactions. The singularity 
ensures the stability of the stable monolayer between surface elevations in Stranski-Krastanov 
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growth. Without wetting interactions the epitaxial film would show Volmer- Weber growth, i.e., 
growth would occur in separated islands not connected by a wetting layer. Note that we here add 
the last term to eq. [4] where fx is the chemical potential. It does not affect the evolution equation 
when M nc = 0. Related (in part non-local) equations are employed in Refs. [|2l l42l - l46l . 

If vapour deposition is used to deposit the epitaxial layer it is to expect that equation ([T]) with 
([4]) and M nc > well describes the process. We are, however, not aware of such an approach in the 
literature. Actually, for large chemical potential fj, as compared to the other terms in F (Eq. ([4])) 
one can even use the system to describe the evolution under vertical deposition of material. The 
constant M nc fi does then corresponds to the constant deposition rate in other models (V in ||38ll ). 

Although the overall form of the equation is identical for the various problems introduced 
above, the specific physics is very different. However, still one can employ the same set of tech- 
niques to analyse the various models. Normally, one uses (i) a linear stability analysis of homo- 
geneous steady states, i.e. flat films to determine the stability of the system and typical length 
scales that will dominate the short-time evolution in case the homogeneous state is unstable, (ii) 
Depending on the properties of the dispersion relation obtained in the linear stability analysis one 
might be able to analytically study stable and unstable steady state solutions and their stability in 
the weakly non-linear regime. This is, however often not possible, (iii) In the strongly non-linear 
regime steady states and their stability might still be obtained, e.g., using continuation techniques 
P71 I48I. These are readily available for two-dimensional systems that can be expressed as ordi- 
nary differential equations Il49l . Recently, they were also introduced for the full three-dimensional 
problem, in particular for Eq. ([T]) with M nc = and ([3]) [50] . Note, that variational methods are apt 
to obtain the steady states directly from the functional F, but are not suitable to discuss the stabil- 
ity of the steady states as this involves dynamic aspects. Many groups prefer to 'skip' step (iii) and 
rather directly (iv) simulate the evolution equation in time using advanced numerical techniques 
(spectral, pseudo-spectral, or semi-implicite). 

We remark here that other types of continuum description exist for all the mentioned systems. 
Whereas here we focus on evolution equations for surface or interface profiles, another class of 
models describes interface evolution using phase fields [1511 . See, for instance, for dewetting 
and liquid films/drops in general Ref. ifTOll and for epitaxial growth Ref. fl2] |43l [52j|. We do en- 
tirely exclude from our consideration the vast literature on discrete stochastic models that exist for 
dewetting/evaporation processes (e.g., Il53l - [55ll as well as for surface growth (for reviews see, e.g., 
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In the following, we restrict ourselfes to two-dimensional physical situations described by film 
thickness profiles <f>(x, t) that depend on one spatial coordinate only. The drops or quantum-dots 
in 2d will actually refer to liquid ridges or quantum-wires in 3d, respectively. We will perform 
a basic analysis of linear stability and steady states 'in parallel' for a dewetting liquid film on a 
solid substrate (Section [TT]) , a dewetting evaporating/condensing liquid film on a solid substrate 



(Section III), and the epitaxial structuring of a solid film (Section IV). Note that we only review 
steps (i) and (iii) of the above introduced scheme that sketches a more complete analysis of the 
system behaviour. The next step would be to use advanced numerical techniques to simulate the 
evolution of the films in time for two- and three-dimensional physical settings. It involves a rather 
large number of techniques and groups and we would like to refer the reader to the individual 
publications cited in the respective sections below. 

Before we start we would like to point out the relevance of stable and unstable steady state 
solutions for systems that evolve in time. Most steady state solutions are either not linearly stable 
or do not correspond to the global energetical minimum that the system will finally approach. Such 
solutions are normally not well appreciated in the literature as they are not present 'in equilibrium' . 
They are, however, often present for a long time in the course of the time evolution and due to their 
character as saddles in function space they do often 'structure' the evolution towards equilibrium. 
Their stability properties (i.e., growth and relaxation rates) determine time-scales for important 
steps of the dynamics, and as important transients in experiments of finite duration they might 
actually even get 'frozen in' or 'dried in' as, e.g., in the dewetting of thin polymer films Il58ll59ll or 
suspensions [|60ll6T| , respectively. In this connection, coarsening, is a particularly interesting issue 
as in its course the system 'passes through' an infinite number of steady state solutions that are 
stable when taking their typical size as reference size, but are unstable with respect to modes on 
larger lengthscales, i.e. with respect to coarsening. One could say the individual solutions do first 
'attract a time-evolution' and then 'expel it' along the single unstable direction (corresponding to 
coarsening). 



II. DEWETTING 

A dewetting film of non- volatile liquid on a solid substrate is modelled by Eqs. ([T]) and ([3]) with 
M nc = [HHHE]. ^ ne l° ca l ener gy f{4>) corresponds to a disjoining pressure 11(0) = — <9^/(0) 
ASHE!- Various functional forms are used for /(</>). The particular choice is, however, not very 
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relevant for the qualitative behaviour of the system. The latter only depends on the number and 
relative depth/height of the extrema of We choose here /(</>) = — n{l / (2(fi 2 ) — 6 3 /(50 5 )] 

as derived by Pismen from a modified Lennard- Jones potential with hard-core repulsion [|62l - |64l . 
The resulting /(</>) has only one minimum at a finite precursO r = b. The parameter k corresponds 
to a typical energy scale. With the chosen signs the first term corresponds to a destabilizing long- 
range van der Waals interaction (i.e., k is proportional to a Hamacker constant) whereas the second 
one represents a short-range stabilizing interaction. In consequence, the model may describe drops 
of a partially wetting liquid in co-existence with a precursor film of thickness b. In the following 
we scale <p by the precursor film thickness, i.e., we fix b = 1. 

If there existed a second minimum at larger finite thickness it would correspond to a criti- 
cal height for the transition between spherical-cap-like drops and pancake-like drops (e.g., the 
capillary length when gravity is included) (6l (65J. For a selection of other pressure terms see, 

e.g. EEnaisaisBisra. 

Inspecting Eq. ([T]) with M nc = one notes that any flat film (thickness <\> = O ) corresponds 
to a steady state solution of the system. However, those films might not be stable: We linearize 
the system about the flat film employing harmonic modes, i.e., <fi(x,t) = O + e exp(/3t + ikx) 
where e <C 1 is the smallness parameter, (3 the growth rate of the harmonic mode of wavenumber 
k. Entering this ansatz in Eq. Q with ([3]) gives the dispersion relation 

P(k) = -M° cl k 2 {k 2 -k 2 c ), (5) 

where M c ° = M c (0 = <p ). The critical wavenumber is given by k c = ■>/— d</></>/| /7- For 
9^/1^=00 < me fil m ls linearly unstable for < k < k c . The most dangerous instability 
mode, i.e., fastest growing mode has /c max = k c /y/2 and /3 max = M c °7^/4. The onset of the 
instability occurs at 0^/]^=^ = with &° nset = 0, i.e., it is a long-wave instability. Note that 
mass conservation implies (3(k = 0) = 0. Examples of dispersion relations above and below the 
instability treshhold are given in Fig. [T] 

Steady film thickness profiles are obtained by solving Eq. ([T]) with d t 4> = 0. Note that in the 
present setting this corresponds to 5F/5<j) = because the first integration constant (when inte- 
grating Eq. [TJ is zero as we look at systems without large-scale meanflow. The whole solution 
structure might be mapped out in parameter space employing continuation techniques using ana- 
lytically or numerically known solutions to start the continuation [fTTl |47] - I49l . Many parameters 
might be used as main continuation parameter, e.g., period or domain size, film thickness, or any 
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wavenumber k 



FIG. 1: Dispersion relations for the instability of a flat liquid film w.r.t. surface modulations resulting in 
dewetting and the evolution of patterns of droplets. Shown is a stable and an unstable case in the generic 
form P = —k 2 (k 2 — k 2 ), i.e., the growth rate [3 is scaled by M c ° 7. The wave number is given in units of k c . 




period period 



FIG. 2: Shown are characteristics of families of steady one-dimensional droplet structures arising in dewet- 
ting of a non- volatile liquid using a simple disjoining pressure. The left panel gives the L2-norm \\5(j)\\ 
and the right one the energy F (Eq. ([3])) per length. The legends give the corresponding mean film thick- 
nesses. The families are obtained using continuation techniques. Note that the solutions are unstable w.r.t. 
coarsening. 
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FIG. 3: Shown are examples of film thickness profiles for the various branches presented in Fig. [2] Mean 
thicknesses (f> are (a) 1.27, (b) 1.5, (c) 2.5, and (d) 5.0. The system size is L = 60 in all cases. Given are 
profiles from the following branches: n = 1 (solid lines), n = 2 (dot-dashed lines), n = 3 (dashed line), 
and n = 4 (dotted line). The thin line in (d) corresponds to a nucleation solution. For given L only the 
n = 1 profiles are stable w.r.t. coarsening. 

parameter of the energy. Here we only use system size. For a more detailed explanation of the 
continuation procedure for thin film equations see the Appendices of 11711 1721 . 

For a linearly unstable film a one parameter family of periodic film profiles bifurcates from the 
flat film at wavenumber k c as illustrated in Fig. |2j At the bifurcation point a steady harmonically 
modulated surface profile exists with a period of 2n/k c and an infinitesimally small amplitude. 
Using period as main continuation parameter the whole family of periodic solutions can be ob- 
tained. A selection of resulting film thickness profiles is given in Fig. [3} We characterize solutions 
by their norm 1 1 50 1 1 = [(l/L) J^((j)(x) — o ) 2 cfe] 1//2 and energy per length E = (1/L) J F[(p]dx. 
Note that mass conservation advises us to directly compare only profiles of identical mean film 
thickness. This is ensured by /i acting as a Lagrange multiplier. In consequence, fx changes along 
the individual branches in such a way that the mean of 4>(x) is kept fixed. 

The resulting branches might bifurcate from the trivial solution supercritically, i.e., they bifur- 
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cate towards the region of the linearly unstable flat film (forward towards larger periods; see, e.g., 
in Fig. [2] curve for O — 1-5). Or the branch bifurcates subcritically, i.e., it emerges towards the 
region of the linearly stable flat film (backward towards smaller periods; see, e.g., in Fig.[2]curves 
for 0o = 1-27, 2.5, or 5.0). The location of the border between sub- and supercritical behaviour 
can be determined employing a weakly nonlinear analysis and can be expressed as an algebraic 
condition for the 2nd, 3rd and 4th derivative of f(<p) 11731 . 

In the case of the subcritical bifurcation, the subcritical branch (i.e., between the bifurcation 
and the saddle-node bifurcation where it 'folds back' towards larger periods) consists of unstable, 
nucleation solutions that aquire an importance for the rupture process of the film if the system is 
noisy or 'dirty': Note that for the chosen potential /(</>) there exist no metastable films. However, 
one may distinguish two thickness ranges within the linearly unstable range: (i) the linearly un- 
stable modes are fast and will dominate the time evolution even in the presence of defects (finite 
size perturbations); (ii) the linear modes are slow, and defects - if there are any present - will 
dominate the time evolution. The distinction is related to the existence of the subcritical branch 
of unstable solutions. They correspond to nucleation or threshold solutions as they have to be 
overcome to break the film into drops smaller than the critical wavelength of the linear instability. 
As they are saddles in function space they can 'organize' the evolution of defect-ridden thin films 
by offering a fast track to film rupture. This allows to determine a typical time for nucleation 
events even inside the linear unstable regime and finally to distinguish the nucleation-dominated 
and instability-dominated behaviour of linearly unstable thin films [[241 1741 . A detailed account 
and comparison to the results of 11731 is found in Ref. 117611 . Note that similar ideas have since 
been applied to the break-up of liquid ribbons [1771 . Recently we also performed a more detailed 
analysis for a three-dimensional systems [[5011781 . 

For unstable flat films of thicknesses <p the respective branches of solutions bifurcating at 
L c = 2n/k c shown in Fig. [2] are only the first of an respective infinite number of primary solution 
branches. These bifurcate at domain sizes L cn = 2irn/k c , n = 1, 2, . . .. The branch bifurcating at 
L cn consist of the n — 1 branch 'stretched in U by a factor n. The actual thickness profiles of the n 
branch consist of n identical drops. This representation of the periodic solutions in dependence of 
domain size might seem pointless for the present problem. The reason is that the different branches 
are entirely decoupled. As the situation changes strongly when either looking at other energy 



functionals (see below Section IV), or when breaking the reflection or translational symmetry of 
the system (by including a lateral driving force [[34l|65l|79l, su bstrate heterogeneity [17111 . or lateral 
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boundary conditions [80]) it is, however, useful to introduce the various branches here. 

The monotonous decrease of energy (on the low energy branch) with system size indicates that 
the system tends to coarse towards larger and larger structure size. The steady solutions on the 
'stable' branches (the branches of drop solutions that continue towards infinite period) are stable 
when looked at in a domain of the size of their period, but unstable on larger domains, i.e., they 
are saddles in function space that form the 'stepping stones' of the coarsening process: They first 
'attract the time-evolution' and then 'expel it' along the only unstable direction (corresponding 
to coarsening). The coarsening behaviour for a dewetting film is discussed in more detail in 
Refs. [|8T| - l84ll . Related results for the Cahn-Hilliard equation are discussed, e.g., in Refs. (31 
MM- 

Note that results on the nonlinear stability of flat films and the related steady state solutions on 
branches not connected to the trivial flat film solutions are not discussed here (but see lfTTll34l f761). 
For discussions of the evolution in time of dewetting films in two- and three-dimensional settings 
we refer the reader to Refs. |0^1El|33lESI5DlEZ5ll72l|S7l|5g|. Note finally that continuation 
may not only be used to obtain families of steady or stationary profiles [[65117011721 173118911 , but as 
well to track their stability and bifurcations Il34ll90ti92ll . 

III. EVAPORATION 

When including evaporation in the presently studied framework one uses Eq. ([T|) with the en- 
ergy functional ([3]) and M nc ^ 0. One notes that unlike the case without evaporation flat films of 
arbitrary thickness = </>o do not correspond to a steady state solution of the system any more. 
However, steady flat film solutions still exist when wettability and evaporation balance. For a 
given chemical potential the corresponding steady thicknesses are given by 

dtp f = n, (6) 

i.e., for the disjoining pressure used here O = [/c(l ± (1 — 46 3 y u/fi;)) 1//2 /2 / u] 1 / 3 . For condensation 
(/i > 0) two such equilibria exist if /i < n/Ab 3 whereas for evaporation (/i < 0) only one thickness 
exists. Note that this might be different for qualitatively different disjoining pressures like the one 
used in [fT3ll . As above the steady films might be unstable. 

We linearize Eq. ([T]) with ([3]) about the flat film of thickness O given by Eq. ([6]) employing 
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FIG. 4: Dispersion relations for the instability of a fiat volatile liquid film w.r.t. surface modulations resulting 
in the evolution of a pattern of droplets. Shown are three qualitatively different cases in the generic form 
= -(k 2 + M° c /M c °) (k 2 - k 2 ), i.e., we fix 7 M C ° = 1. We furthermore choose M n ° c /M c ° = 1/4. 

harmonic modes as above to obtain the dispersion relation 

0(k) = ~M Q cl {k 2 + ^j (k 2 -k 2 ) (7) 

where M° c = M nc (0 = O ), M c ° = M c (0 = <j> ), and k c = y / -d (j)< t > f\ 4)=<t>0 />y as without evapora- 
tion. As the expression in the first parenthesis is always positive the film is linearly unstable for 
9^/1^=00 < in the wave number range < k < k c . The most dangerous instability mode occurs 
for k mSLX = ^/(k 2 - M n ° c /M c °)/2 with /3 max = M c ° 7 (M n ° c /M c ° + k 2 ) 2 /4 if the expression under 
the square root is positive. Otherwise the mode with k = grows fastest (f3(k = 0) = M® c k 2 ). 

The onset of the instability occurs at 9^,/|^=^ = with k° nsct = 0, i.e., it is a long-wave 
instability. Qualitatively different examples of dispersion relations above and below the instability 
treshold are given in Fig. |4j 

As above steady thickness profiles are obtained by setting d t (f) = in Eq. ([T]). For a linearly 
unstable film a one parameter family of profiles bifurcates from the fiat film at k c as illustrated 
in the four panels of Fig. |5j In contrast to the case without evaporation/condensation here the 
mean film thickness changes along the branches whereas the imposed chemical potential remains 
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FIG. 5: Shown are characteristics of families of steady one-dimensional droplet structures arising in dewet- 
ting of a volatile liquid using a simple disjoining pressure. From top left to bottom right the panels give 
the Z/2-norm, the energy F (Eq. (|3]>) per length, the mean film thickness <p, and the droplet amplitude. The 
legends give the corresponding mean film thicknesses for the flat film (po at the identical chmicalpoten- 
tial jjL. The families are obtained using continuation techniques. Note that the solutions are unstable w.r.t. 
coarsening. 

constant. A selection of film thickness profiles is given in Fig. [6] 

Most notably the resulting branches always bifurcate supercritically from the trivial solution. 
The results furthermore indicate that there is no metastable film thickness range at all, implying 
that nucleation does not play any role for a condensing/evaporating film described by the present 
disjoining pressure. This, however, needs a deeper analysis in the future including a comparison of 
the behaviour for different pressure terms. Remarkably, the amplitude (droplet height) approaches 
a limiting value for a domain size of about 3L c /2. As the energy monotonously decreases with 
domain size we expect these solutions to be unstable w.r.t. coarsening as the droplet pattern in the 
non-volatile case. Note finally, that everything discussed above for non-volatile films regarding 
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FIG. 6: Shown are examples of film thickness profiles for the various branches presented in Fig. [5] Mean 
thicknesses (f> are (a) 1.27, (b) 1.5, (c) 2.5, and (d) 3.0. The system size is L = 60 in all cases. Given are 
profiles from the following branches: n = 1 (solid lines), n = 2 (dot-dashed lines), n = 3 (dashed line), 
and n = 4 (dotted line). 

the infinite number of primary solution branches as well applies to the branches shown in Fig. [5] 

After the short analysis of the the volatile and non- volatile thin liquid film system we next focus 
on a thin film equation that describes the evolution of a solid film. 



IV. EPITAXIAL GROWTH 



The final example is a thin film model for the epitaxial evolution of surfaces of crystalline solids 
for a fixed ammount of deposited material that is illustrated here employing a simplified 'glued 
wetting-layer model' [|40l . The evolution equation ([T|) describes such a system when combined 
with the energy functional Q and M nc = 0. As in the case of the non- volatile liquid film material 
is conserved and any fiat film of thickness <p — 0o corresponds to a steady state solution of the 
system. 

As those trivial steady states might be unstable we perform a linear stability analysis along the 
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FIG. 7: Dispersion relations for the instability of a flat epitaxial film w.r.t. surface modulations resulting 
in the growth of quantum dots. Used is the simplest model proposed in Ref. [ 40 ] . The legend gives the 
dimensionless value of d^f at Starting with the lowest one the three bold curves correspond to the cases 
studied in Figs. 4(a), (b) and (c) of [40]. 

lines of the previous sections and obtain the dispersion relation 

/3(k) = -M> k 2 (k 2 - k 2 cl ) (k 2 - k 2 2 ) (8) 

with 

/.-. 



\ 



I (9 ) 



a 2 



"el/c2 

The film is linearly unstable for d^f]^^ < a 2 /Av in the wave number range k c \ < k < k c2 . 
The fastest growing mode has the wavenumber k max = (a /Sv) 1 / 2 \J I + — "&vd^f\ /o^. 
For values of d^f]^^ < one of the two critical wavenumbers becomes imaginary, i.e., the film 
is then unstable for < k < k c2 . Note that the onset of the instability occurs at d^f]^^ = a 2 /Av 
with k° Dset = yfaj2v, i.e., it is a short-wave instability. Qualitatively different examples of actual 
dispersion curves above and below the instability treshold are given in Fig. [7j 
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FIG. 8: Shown are characteristics of families of steady one-dimensional surface structures (quantum dots) 
arising in epitaxial growth using the simplest model proposed in Ref. [40]. The panels in the left column 
give the L2-norm and the ones in the right column the energy F (Eq. Q) per length. Parameter values 
correspond to (top row) Fig. 4(a), (middle row) Fig. 4(b), and (bottom row) Fig. 4(c) of [40]. From top to 
bottom the non-dimensional wetting interaction increases leading to more intricate behaviour. The families 
are obtained using continuation techniques. 
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FIG. 9: Shown are examples of film thickness profiles (<p — phio) for the various branches presented in 
Fig. [8] The parameters correspond to (a) the top row, (b) the middle row, and (c) the botom row of Fig. [8] 
The chosen system size is L = 17 in all cases. Given are profiles from the n = 1 + branch (solid lines), the 
n = 2 + branch (dot-dashed lines), the n = 1~ branch (dashed line in (b)), and the first side branch of the 
n = 2 + branch (dotted line in (c)) as defined in Fig. [8j Mean film thicknesses is 0o = 1.75. 

As in the other two cases steady thickness profiles are obtained by setting d t <p = in Eq. ([T]). 
Here they do not represent profiles of liquid droplets, but profiles of solid quantum- or nanodots. 
Due to the higher order terms in the energy ([4]) as compared to the functional ([3]) studied up to now, 
the dispersion curves (Fig. [7]) are qualitatively quite different from the ones for dewetting (Fig. [T]). 
In consequence we expect a rather different picture for the steady state solutions as well. At the 
parameter values allowing for two critical wavenumbers k cl and k c2 , i.e. for d^f]^^ > 0, one 
finds two bifurcations from the trivial solutions at domain sizes 2tt /k c \ and 27r//c c2 , respectively. 
Near onset one expects the two emerging branches to interconnect, as at onset the solution space 
should change only locally. Further above onset the behaviour might, however, be dramatically 
different. The two branches bifurcating at 2k jk cl and 27r//c c2 are the first of two infinite series of 
primary solution branches that bifurcate at respective domain sizes L cln = 2nn/k cl and L c2n = 

2%n/k C 2 with n = 1,2, This consideration is of importance here, as the various primary 

branches may actually couple. 

We illustrate the solution structure in Fig. [8] using parameter values as in Fig. 4 of Ref. iPlOll . 
Shown are the norm (left column) and energy (right column) of the solutions on the various 
branches for three (positive) values of the strength of the wetting interaction. A selection of film 
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thickness profiles is given in Fig. [9]. As the strength of the wetting interaction decreases the system 
gets deeper into the unstable regime and the behaviour becomes more intricate. In the upper row 
of Fig. [8] we are slightly above onset and the two n = 1 branches that emerge at the two zero 
crossings of the dispersion relation actually connect as expected. They are well separated from 
the branches with n > 1. Single-quantum-dot (n = 1) solutions only exist for a small range of 
periods L cl < L < L c2 , i.e., the system is not able to undergo any coarsening. This corresponds 
well to the results obtained in ll40l using weakly nonlinear analysis and integration in time. In the 
following we will call the branch emerging at iL cl [iL c2 ] the left [right] n = i branch. 

The picture changes, however, further above onset, i.e., when decreasing the wetting interac- 
tons (middle row of Fig. [8]). The two n = 1 branches do not interconnect any more. Actually, the 
left n = 1 branch continues towards L = oo without any side branches. The left n = 2 branch may 
change its stability with respect to coarsening with increasing domain size. Every time it changes 
stability a secondary branch bifurcates in a period doubling bifurcation rather similar to the sce- 
nario described for driven liquid films on inlined plates Il34l . The right n = 1 branch connects 
via one of the branching points to the left n = 2 branch. The behaviour becomes increasingly 
involved for higher branch numbers. 

Decreasing the wetting interaction even more brings the system further above threshhold (bot- 
tom row of Fig. [8]). As k c \ becomes quite small the right n = 1 branch bifurcates at very large 
domain sizes. However, the first secondary branch of the left n = 2 branch now split up into 
various pieces that connect different left branches via secondary bifurcations. 

We will at present not go deeper into the analysis. In particular no proper stability analysis for 
the steady states is done here and as above for the dewetting system time evolution is not touched. 
Simulations can be found, for instance, in Refs. ll2l l4{Hl42ll44l - l46ll . Note that the results on stability 
sketched in the previous paragraph are infered from the structure of the bifurcation diagrams alone. 
They need to be scrutinized in detail as they are very important for the coarsening behaviour. The 
non-monotonous dependence of energy on system size in Fig. [8] indicates that coarsening will 
depend on structure length in a non-trivial way. The long-time corsening for epitaxial growth is 
discussed, e.g., in Refs. ll93H95i 

The results presented for the glued wetting-layer model indicate that our knowledge about 
epitaxial growth would benefit from a more detailed analysis that maps out the steady solutions 
and their stability systematically for the various models proposed in the literature. Numerical 
techniques are available to do this not only for two-dimensional but as well for three-dimensional 
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systems [50J. 

V. CONCLUSION 

In the present contribution we have applied a basic mathematical analysis to three different 
physical systems involving solid and liquid films at solid surfaces that may undergo a structuring 
process by dewetting, evaporation/condensation or epitaxial growth, respectively. The aim has 
been to highlight similarities and differences of the three systems based on the observation that 
all of them can be described using models of similar form, i.e., a time evolution equation for a 
purely dissipative system without any inertia based on a gradient dynamics that is characterized by 
mobility functions and an underlying energy functional. Thereby the dynamics might have a non- 
conserved and/or conserved part. Equations, like the well known Cahn-Hilliard (purely conserved) 
and Allen-Cahn (purely non-conserved) equations represent limiting cases (with respect to the 
chosen dynamics) of this general evolution equation. They are normally used with a particularly 
simple form of the energy functional as well. 

The three physical systems discussed here all fit into the general framework when chosing 
different mobilities and/or energy functionals. Beside the mathematical similarity the systems do 
as well model similar physical phenomena as the respective liquid and solid films structure under 
the influence of their effective interaction with the substrate. However, although most readers will 
agree that the epitaxial growth of quantum dots is a subject of surface science less of them might 
do so in the case of a dewetting liquid film. We hope that our review contributes to the development 
of a more unified view onto these processes of self-organisation at interfaces. 

In particular, we have used two basic steps of mathematical analysis, namely the study of the 
linear stability of homogeneous steady states, i.e., flat films, and the mapping of non-trivial steady 
states, i.e., drop/hole (quantum dot) in dependence of system size for various values of interaction 
constants and/or mean film thickness. Our aim has been to illustrate that the underlying solution 
structure might be very complex as in the case of epitaxial growth but can be understood better 
when comparing to the much simpler results for the dewetting liquid film. We have furthermore 
shown that the continuation techniques employed can shed some light on this structure in a more 
convenient way than time-stepping methods. One might further argue that understanding the solu- 
tion structure of the quantum-dot system might as well allow to predict pathways of evolution in 
time for an epitaxial layer as the behaviour will ultimately be determined by the steady solutions 
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and their stability. This procedure was already followed for thin liquid films f7Tll74ll96ll . 

The usage of a general formulation like Eq. ([TJ) does not only relate seemingly not related phys- 
ical systems mathematically. It does as well allow to discuss model extensions in a more unified 
way. Taking epitaxial growth as example, the general form ([[]) would propose that in the non- 
conserved case, i.e., when material is deposited from the gas phase continuously the non-conserved 
part of the equation should be — M nc 5F/5(p with F given by ([4]). The constant depostion often 
used in the literature would then only refer to the limit of a rather large chemical potential. Oblique 
deposition of material Il38ll might as well be modelled incorporating beside the vertical influx as 
well lateral driving contributions. 

The here reviewed formulation can be extended in a straightforward way to the case of N 
coupled fields 4> = (0i, 2 , . . . , <p^) following a gradient dynamics that is governed by an energy 
functional F[cpi, (f> 2 , . . . , 0at]- The evolution equation is 

5F 



d t cf> = V- 



SF 



M ra ■ ^ (10) 



with the M c (0) > and M ne (0) > being symmetric positive-definite mobility matrices for 
the conserved and non-conserved part of the dynamics, respectively. They are formed by N x N 
mobility functions, respectively. Note that 5F/5<f> corresponds to a vector of dimension N. A 
typical example is the evolution of a dewetting two-layer system where the <pi and 2 correspond 
to the two film thicknesses, respectively. The formulation as a special case of ( flQ] ) is derived 
from the Navier-Stokes equations in the two liquid layers and appropriate boundary conditions for 
spatially two- and three-dimensional settings in Refs. [|96] - l99l . The formulation in the form of 
Eq. ( [10] ) has the advantage that one can easily check the consistencey of the mobility functions 
after a normally rather involved derivation. Furthermore one can draw on rather general results, 
e.g., for the linear stability of flat film solutions 11971 . Note, that such a formulation will apply to 
any multi-layer system in a relaxational setting, including films of dielectric liquids in a capacitor 
HI 001 1101H or multilayer pending films under the influence of gravity. We expect it to hold as well 
for the relaxation dynamics of multilayer epitaxial films. 

Let us finally note that we have entirely excluded spatially heterogeneous systems, i.e., systems 
described by equations of form ([TJ where the energy depends explicitly on position. Such systems 
are already treated for dewetting films [[171 17T1 [821 1102J - I104L however, we are not aware of any 
such studies for epitaxial growth. We have as well excluded laterally driven systems described 
by equations like ([TJ that additionally include lateral driving forces. A comparison of sliding 
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droplets on an incline ll34l l88l . a driven (or convective) Cahn-Hilliard system IT791 and epitaxial 
growth under oblique deposition might be rather interesting, for instance, in terms of coarsening 
behaviour [|34l |94l [T05l 110611 . The combination of both effects, i.e., heterogeneous systems with 
lateral driving allows to describe the depinning and stick-slip motion of droplets ll92l I107H . The 
problems has up to now, however, no counterpart in epitaxial growth. 
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